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Abstract 

The development of any super/hypersonic aircraft 
requires the interaction of a wide variety of technical 
disciplines to maximize vehicle performance. For flight 
and engine control system design and development on 
this class of vehicle, realistic mathematical simulation 
models of atmospheric turbulence, including winds and 
the varying thermodynamic properties of the 
atmosphere, are needed. A model which has been 
tentatively selected by a govemment/industry group of 
flight and engine/inlet controls representatives working 
on the High Speed Civil Transport is one based on the 
Kolmogorov spectrum function. This report compares 
the Dryden and Kolmogorov turbulence forms, and 
describes enhancements that add functionality to the 
selected Kolmogorov model. These added features are: 
an altitude variation of the eddy dissipation rate based on 
Dryden data, the mapping of the eddy dissipation rate 
database onto a regular latitude and longitude grid, a 
method to account for flight at large vehicle attitude 
angles, and a procedure for transitioning smoothly 
across turbulence segments. 

Sy mbols 

A„ Fourier series cosine coefficient 

B n Fourier series sine coefficient 

E Kolmogorov energy spectrum 

fi lower bound on frequency range of interest, Hz 

f 2 upper bound on frequency range of interest, Hz 

h altitude, km 

kj lower bound on wavenumber range of interest, 

rad/m 

k 2 upper bound on wavenumber range of interest, 
rad/m 

L integral scale length, m 

Lq Kolmogorov outer scale length, m 

u longitudinal turbulence velocity, m/s 

v transverse turbulence velocity component, m/s 
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V air relative vehicle velocity, excluding 
turbulence components, m/s 
w turbulence velocity component orthogonal to u 
and v, m/s 

x distance traveled along flight path, m 

a general Kolmogorov energy spectrum constant 

a p Kolmogorov pressure energy spectrum 

constant, hPa 2 s 2 /m 2 
y flight path angle, deg 

Ap density perturbation from mean 

AP pressure perturbation from mean 

AT temperature perturbation from mean 

£ eddy dissipation rate, m 2 /s 3 

e D eddy dissipation rate derived from Dryden data, 
m 2 /s 3 

£ 0 eddy dissipation rate used in A,, and B„ 
calculation, m 2 /s 3 
O Dryden energy spectrum 

O v von Karman energy spectrum 

§ body axis roll angle, deg 

0 pitch attitude, deg 

y w heading relative to wind, deg 

Cl wavenumber, rad/m 

a RMS turbulence velocity, m/s 

Ab b r e _ viatiQn$ 

3-D three dimensional 

2-D two dimensional 

lat. latitude 

long. longitude 

HSCT High Speed Civil Transport 

RMS Root Mean Square 

Su bs cri pts 

h horizontal component 

u longitudinal component 

v transverse component 

w component orthogonal to u and v 

0 indicates ffeestream value 

Introduction 

The development of any new aircraft requires the 
interaction of a wide variety of technical disciplines. 
Furthermore, when the vehicle being designed requires 
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the development of new technology in many of these 
disciplines, the level of interaction between the s 
technology teams must be high to maximize efficiency 
in the resulting design. Such is the case with the 
technology development of the High Speed Civil 
Transport (HSCT) under the High Speed Research 
program. 

In the case of the HSCT program it was recognized that 
different disciplines used different atmosphere models. 

In order to establish a common basis for comparison 
and minimize duplication of effort, an ad hoc group, 
including representatives from NASA and industry, has 
formed with the goal of adopting a model of 
atmospheric turbulence. The selected model will be 
used within the engine community as a basis for the 
evaluation of inlet/engine performance and in the flight 
controls area for the evaluation of the aircraft control 
system and handling qualities in piloted and batch 
computer simulations. 

The model which has been proposed by the group is one 
suggested by Bill Tank of Boeing 1 . The model, like the 
von Karman model, is defined in the frequency domain, 
and for which the power has a -5/3 intensity roll off 
with increasing frequency. A key reason for the 
selection of the model is the capability to generate 
thermodynamic perturbations (density, pressure, and 
temperature) as well as wind gusts. A second reason for 
selecting the Tank model is the data used to 
parameterize it. Tank bases the model on data collected 
from a variety of sources where data was taken in flight 
as close as possible to the approximate cruise altitude of 
the HSCT. 

This report reviews the Tank model, discusses 
implementation issues, and presents extensions made to 
the model to enhance its capability. 

Background 

The atmosphere has long been studied to better 
understand and predict its behavior, and for aircraft 
designers and operators, understanding the turbulent 
nature of the atmosphere is an area of concern. Current 
atmospheric theory 2 states that the turbulent energy of 
the atmosphere is distributed as a continuous function 
of wavenumber (spatial frequency) as shown in Figure 
1 . This Figure shows the zonal (east- west) wind, 
meridional (north-south) wind, and potential temperature 
power spectra (note: for clarity the meridional and 
temperature spectra have been shifted one and two 
decades to the right respectively). Notice that at the 
very long (planetary scale) wavelengths the power of 
turbulence decreases as the -3 power of decreasing 


wavelength (increasing frequency ), and from about a 
wavelength of 500 km, the power fall off is roughly as 
the -5/3 power of frequency. Turbulence theory also 
states that the range where turbulence is fully three 
dimensional and isotropic extends from the very small 
scale (< 1cm) out to a scale of a few hundred meters 3 . 

In the range between a few hundred meters and 500 km 
the -5/3 law is still obeyed for turbulence in the near 
horizontal plane, though the turbulence cannot be 
classified as three dimensional (3-D) isotropic. Vertical 
turbulence has been found to obey a power fall off of -3 
to -2 down to scales of 100 meters. It is expected to 
change to a -5/3 slope at even shorter scales 3 . 

In the 3-D region turbulence energy is not so much lost 
due to friction as it is transferred to shorter and shorter 
wavelengths. The rate at which energy is cascaded down 
in scale is called the eddy dissipation rate, £. There is 
also evidence that supports the theory that energy is not 
only cascaded down in scale but also up. That the -5/3 
law holds for scales out to hundreds of kilometers is a 
fortunate result of upscale transfer of two dimensional 
(2-D) turbulence energy. The impact this has on the 
modeling of turbulence for simulations is that models 
developed for 3-D, isotropic turbulence at short scales 
can be extended beyond the 3-D regime into the margins 
of the 2-D without changing the model. 

In order to evaluate the performance of aircraft systems 
and flying qualities an appropriate, relatively simple, 
characterization of atmospheric turbulence is required. It 
is not sufficient to always assume a quiescent 
atmosphere or even that the thermodynamic properties 
of the atmosphere are solely deterministic functions of 
altitude. 


Traditionally the Dryden turbulence model has been used 
whenever atmospheric (wind) turbulence was needed to 
add some sense of realism to a time-based simulation or 
to provide a disturbance input when evaluating systems 
in the frequency domain. The Dryden spectrum for the 
longitudinal, body axis velocity disturbance is given 4 
as, 


0(O) U = o' 


2L„ 


1 


0 ) 


* i+(L u nr 

where Q is the wavenumber, or frequency term, in units 
of radians per meter, a represents the standard deviation 
of the velocity disturbance, and L u is the integral scale 
of the longitudinal turbulence, which will be discussed 
in more detail later. 


The spectra of the lateral and vertical components of 
Dryden body axis turbulence are given in the same 
reference by, 
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with a similar definition of terms. Although the 
Dryden equations are more frequently used, they are not 
the best fit to the observed character of the atmosphere. 
Notice in equations 1 and 2, that the Dryden spectra fall 
off as the -2 power of frequency, not the -5/3 power as 
shown in Figure 1 . The predominant use of the Dryden 
form stems from the ease of implementation within a 
time-based simulation. Such implementation usually 
involves forming a band limited, white noise filter 
equation by taking the spectral square roots of the 
spectrum equations. This is not so simple for an 
irrational -5/3 power fall off without approximating, 
but is relatively simple for the rational -2 power fall in 
equations 1 and 2, The resulting transfer functions 4 are 
given below, 
longitudinal: 
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lateral and vertical: 
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Use of these equations with a band-limited, unit 
variance, white noise input and appropriate filter 
parameters yields Dryden aircraft-scale turbulence. For 
the puipose of classical rigid-body simulation and 
piloted assessment, the Dryden form is sufficient. 

When aircraft flexibility is an issue for vehicle flight 
assessment the frequency range of interest widens to 
include the higher frequency structural modes. At this 
point the difference between the -2 and the -5/3 power 
fall off can not be ignored. 

The spectrum functions that better fit experimentally 
observed turbulence are the von Karman spectra as given 
in the same reference 4 . 
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for the lateral and vertical components. 

Both the Diyden form and the von Karman are 
parameterized by two terms, the standard deviation of 
the turbulence, a, and an integral scale length. The 
standard deviation is equal in magnitude to the square 
root of the area under a power spectrum plot. Note that 
increasing or decreasing a in the equations 1, 2, 5, or 6 
would have the effect of shifting up or down the 
amplitude of the spectrum, without distorting the shape 
of the plot. The Dryden spectra shown in Figure 2 were 
generated with the equation, 




10 

1 + (0.5Q) 2 


(7) 


Referring to eqn. 1, this equation arises when L u = 0.5m 
and c 2 = IOtc m 2 /s 2 . The solid line in the upper plot on 
Figure 2 shows the typical Dryden spectral shape on a 
log-log plot. Note that the log-log plot is deceptive 
when attempting to characterize the area distribution by 
visual inspection. Although most of the area, and thus 
larger contributor to the integral of the spectrum, 
appears to lie under the flat (low frequency) portion of 
the power spectrum, that relation is not correct. As 
decade by decade integration of the spectrum would 
reveal, most of the turbulence power being described by 
the plot, is under the sloping portion of the curve. 


The upper plot is a log-log plot and the lower a semi- 
log of the same data, with the solid lines defined by 
equation 7. The dashed line in Figure 2 represents the 
product of the power (or Dryden amplitude) and the axis 
frequency. Plotting this product as a function frequency 
gives a more realistic (less deceptive) graphical 
representation of the area distribution. 

The knee, or the point where the spectrum has fallen to 
the half power point with increasing wavenumber (or 
frequency), is located at a wavenumber of 1/L. In 
equation 7, L has a value of 0.5. As shown in Figure 
2, the knee is properly located at a wavenumber of 1/L U , 
or 2.0 Rad./m. For proper use, the wavenumbers of 
interest to the designer should be near 1/L U , so that the 
power associated with those wavenumbers is 
appropriate. Arbitrarily increasing the scale length, L, 
has the effect of shifting the knee of the spectrum to the 
left. If the knee moves left and a is held constant, the 
spectrum will be shifted left and upward slightly to keep 
the integrated area constant. This shifting reduces the 
power associated with the higher frequencies. In order to 
both not affect the higher frequency power distribution 
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and move the knee to the left (for example, to include 
more low frequency power), both L and 0 must x 
increase. So it should seem intuitive now that the 
standard deviation and scale length should be related. 
Tables defining L and 0 pairs, exist and give the pairs 
as a function of altitude, such as is shown in Table l 5 . 

Despite being frequently used there is much controversy 
over specification of the scale length. According to 
Etkin, the determination of the scale length at altitude 
remains an unresolved difficulty. This stems from the 
uncertainty of how to determine L, whether from the 
integral of the velocity autocorrelation function or from 
a fit to a spectrum 6 . The matter is further confused by 
the tendency to think of the integral scale as a physical 
characteristic of the turbulence field, and is often 
confused with the ‘outer scale* of turbulence, or the 
longest wavelength that could be expected. For 
example, in Figure 2 the integral scale is the reciprocal 
of the knee frequency, or 0.5, but inspection of the bell 
shaped area distribution indicates that significant 
turbulence intensity exists at longer scales (lower 
frequencies). 

Additionally, the use of the Dryden or von Karman 
models provide turbulent winds, but not thermodynamic 
variations. The models do not usually account for the 
global variation of turbulence. In time-based 
simulations, where the Dryden model is used, the power 
fall off is not the observed -5/3 value, though some 
may argue that the Dryden model gives a better fit to 
the vertical gusts due their observed -3 to -2 character at 
scales above 100m as mentioned above. 

The Tank Turbulence Model 
In the light of these issues with using the Dryden and 
von Karman models, it is interesting to note that 
atmospheric scientists use a different model. Where the 
aerodynamicist uses the Dryden or von Karman models 
for characterizations of the turbulent atmosphere, 
atmospheric scientists tend to use the Kolmogorov 
energy spectrum 7 . The Kolmogorov spectrum function 
for longitudinal turbulence has the form, 

2 -5 

E(k) = oce 3 k 3 (8) 


The turbulence model suggested by Tank 1 for 
application to the HSCT, hereafter referred to as the 
Tank model, uses this Kolmogorov form for defining 
wind and thermodynamic turbulence. According to the 
Tank model the temperature spectrum has the same 
form as the longitudinal, the only difference being the 
constant a, whose units vary depending on the 
spectrum being defined. A description of the similarity 
of the temperature and velocity spectra is available 9 . 

The spectrum of atmospheric pressure perturbations 
associated with turbulence is given by Tank to be 

E(k) = a p (P/T 0 ) 2 e 273 k* 5 ' 3 (9) 

where P () is the standard atmospheric pressure in hPa and 
T () is the standard Kelvin temperature. 

The Tank model defines the temperature and pressure 
equations, leaving the the density perturbation to be 
determined from the equation of state. If, 

(P 0 + AP) = (p 0 +Ap) R (T 0 + AT) (10) 


where p 0 is the standard atmospheric density, R is the 
gas constant, and the A symbol indicates a perturbation 
of an atmospheric quantity (temperature, T, or pressure, 
P), then the density perturbation can be solved for 


(Po+AP) Po_ 

(To + AT)R RTo 


(ID 


Given these equations and the appropriate a value for 
the quantity being simulated, the Tank model reduces to 
one parameter, e. This contrasts with conventional 
Dryden and von Karman which have two parameters (a 
and L). The eddy dissipation rate, £, characterizes the 
turbulence for whatever range of frequencies is of 
interest. To account for the global non-homogeneous 
nature of turbulence. Tank has defined e, at the 
approximate HSCT cruise altitude, as shown in Figure 
3. Inspection of the figure reveals that Tank has defined 
£ between 80 degrees north and 40 degrees south latitude 
over the oceans around the world. 


The eddy dissipation rate can be linked to more familiar 
terms by integrating both sides of eqn. 8 from kj to k 2 
and recognizing that the integral of the left side is equal 
to 0 2 . This yields the following equation, 


where k is the spatial wavenumber, £ is the eddy 
dissipation rate, and a is a specified constant. In log 
space eqn. 8 is the equation for a straight line with a 
-5/3 slope. Notice that it matches the -5/3 character 
observed in Figure 1, and is similar to the von Karman 
form in the limit as wavenumber approaches infinity. 

In the 3-D range the transverse spectra are 4/3 times the 
longitudinal for an incompressible fluid as Kolmogorov 
showed 8 . 
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Here a 2 represents the area under the spectrum defined in 
eqn. 8 and is the variance of the turbulence within the 
range k ( to k 2 . Integrating from k, to infinite 
wavenumber equation 12 reduces to 

<7 3 

£ = 3 03) 

(1.5a)5Lo 

where Lq is 1/ k t and represents the longest wavelength 
in the range of interest and is referred to as the outer 
scale length. This outer scale should not be confused 
with the integral scale length. 

Tank Model Extensions 

To implement and add functionality to a time-based 
simulation employing the Tank model, several features 
have been developed by the author in a series of 
MATLAB scripts. A variation of the eddy dissipation 
rate with altitude based on Dryden data has been 
calculated. The e database has been modified. Methods 
to more properly account for flight at large attitude 
angles have been implemented. Also, a procedure for 
transitioning smoothly from one turbulence segment to 
the next has been developed. 

Altitud e V a riati on 

The model proposed by Tank defined the eddy 
dissipation rate only at the estimated cruise altitude 
(—60 kft) as shown in Figure 3. To use the model in a 
broader range of the flight envelope a method for 
changing e based on altitude was developed using eqn. 

13 and data from reference 3. 

Table l 5 defines a set of standard deviations and integral 
scale lengths for light, moderate, and severe turbulence 
as a function of altitude for use in a Dryden turbulence 
model. Using this a and L information for horizontal 
turbulence, with a (in eqn. 13) set equal to 0.15 as 
specified by Tank for longitudinal turbulence, eddy 
dissipation rates were estimated from the Dryden 
parameters. These estimated eddy dissipation rates are 
shown in Table 2. 


observation is substantiated by the overlap of the power 
spectra at the low frequencies. Also note that the Tank 
model produces the higher amplitude high frequency 
turbulence, this is due to the difference in power fall off 
between the two models. As should be expected by 
looking at the power spectra the standard deviation of 
the Tank time history is higher than the Dryden (.2699 
vs. .2158 m/s). 

Accepting the similarity, the light turbulence £ data was 
normalized to unity at 18 km and is shown in the last 
column of Table 2. The resulting normalized variation 
with altitude is the altitude variation that has been added 
through eqn. 14, 

e(lat., long., h) = e(lat., long.) * e D (h) (14) 

where h is altitude in km. 

Da ta b ase jGnd difl g 

The £ data given by Tank 1 was only presented in the 
form of the isoline plot shown as Figure 3. For 
implementation within a computer simulation, the data 
was mapped into a regular latitude and longitude grid as 
shown in Table 3. The missing data over the 
continents was filled in with a constant value of 
.000021 corresponding to the global average £ value of 
.00002 plus an offset (the insignificant offset was used 
for bookkeeping purposes to identify the value as one 
not originally defined by Tank). 

La r ge An gl$ Flig ht 

To explain the remainder of the additions to the Tank 
model first requires the explanation of the time-based 
implementation method. 

With the spectrum defined in eqn. 8, any of a number of 
methods could be used to evolve a time simulation. 
Tank suggests a method described by Fairall, White, and 
Thompson 3 . If the spectrum is discretized in M 
frequency bands Ak wide between kj and k 2 , the strip 
area can be computed and equated as 

area = 0.5 * (A„ 2 + B„ 2 ) (15) 


Noting that Tank defines the lower stratosphere global 
mean epsilon to 0.00002 m 2 /s 3 and observing that at an 
altitude of 18 km the light, Dryden turbulence data 
produced a nearly identical epsilon value (shaded in 
Table 2), it was accepted that the iighf turbulence data 
best fit Tank’s data. To support that claim, two time 
histories were evolved, one by the Tank method and the 
other with the Dryden filter. Sample time histories can 
be compared in Figure 4 and the corresponding power 
spectra in Figure 5. Notice that both models generate 
nearly the same low frequency character and that this 


where Ak is defined as numerically equal to kj (the 
lowest frequency of interest), and n ranges from 1 to M. 


After solving for A n and B n , with a uniform random 
split between A n and B n , a Fourier series can be used to 
represent a turbulence field as follows 


M 

u(x)= I A n cos 

n=l 


'27txn' 

+ B n sin 

r 2^xn^ 

< Lq j 


< > 


(16) 


Here L,, is the segment length 1/k,. When the sequence 
above is evolved a one dimensional turbulence field in 
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space is defined. Using Taylor’s frozen field hypothesis 
the space field can be converted to a time-based sequence 
with the simple transformation 

x = V * t (17) 

Incorporating this and exchanging Ak for 1/ Ly , eqn. 16 
becomes 
M 

u(t)= 1 A n cos(2Mktn) + B n sin(2flAktn) (18) 

n=l 

If evolved for a time longer than it takes to travel a 
distance Lo, as defined by eqn. 17, the field will begin to 
repeat. To avoid this a new set of A n and B n coefficients 
must be computed and used each time a distance Lq is 
traversed. 

Given this method of generating a turbulence field, 
several comments can be made. The frequency range of 
interest is defined in terms of spatial waves (if the 
vehicle velocity changes, the effective frequency 
changes). The turbulence field is evolved in 
precomputed segments. The vehicle state does not 
continuously update the model. The vehicle orientation 
has no effect on the sensed winds. 

Dissipation Rate Scaling 

To account for the fact that e may change significantly 
during a segment that maybe as long as 100 km, £ 
scaling has been implemented. After computing the 
turbulence field with a fixed eddy dissipation rate, £(„ 
the instantaneous value of e can be used to form a 
scaling term, 

scale factor = ( e/Eq ) 1/3 (19) 

Equation 19 was derived through algebraic manipulation 
of equation 15 (with the strip area equal to the integral 
of equation 8) and observing the effect of E on A n and 
B n . This scaling factor is applied directly to all the 
turbulence quantities. Figure 6 contains two traces of 
temperature as a function of time. The upper trace 
shows dissipation rate scaling as £ is arbitrarily 
increased by a factor of 4 over the 250 second run. The 
lower trace, which has been shifted down 2 degrees for 
clarity, is the same sequence without scaling. Notice 
how the magnitude of the perturbations increases in the 
upper trace and not in the lower. 

Coordinate Sy stem Rotation 

As in the case of Dryden turbulence, the Tank model as 
defined produces body axis turbulence for horizontally 
flying vehicles. This restriction is the result of the 
turbulence filters being derived for a turbulence field 
aligned along and normal to the vehicle velocity vector. 
Also, for horizontal flight, the vehicle altitude is nearly 


constant so the filter parameters, which are usually 
functions of altitude, do not change. To allow 
simulation of flight at large vehicle attitude angles, 
where there may be significant misalignment of the 
velocity vector and x-body axis, or large roll angle, a 
transformation has been defined. 

To develop the axis transformation equations, the axis 
systems had to be defined. For ease of implementation 
it was desired that the output of the turbulence model 
should be body axis turbulence. It was also known that 
the longitudinal turbulence equation was specifically 
derived to generate turbulence along the velocity vector. 
The two remaining axes are only required to be normal 
to the longitudinal (and each other). For simplicity, in 
this implementation, one of these axes always lies in 
the local horizontal. The third axis points as required to 
remain normal to the local horizontal and the velocity 
vector. This choice has the benefit of aligning with the 
gravity and body axes in the case of standard horizontal 
flight. 

To compute the transformation the heading ‘off the 
wind’ 0F W ), pitch attitude (0), flight path angle (y), and 
roll angle (<J>) are required. Rotating from these 
turbulence axes to body axes involves first a rotation 
through *F W about the local vertical axis. Second, is an 
upward rotation about the one axis lying in the local 
horizontal plane that is normal to the velocity vector, 
through the angle 0- y. The final rotation is about the 
x-body axis through the roll angle <|>. 

The turbulence to body axis transformation matrix can 
be written as the reverse order multiplication of the 
individual vector rotations as follows 

t BT ■ 


’l 

0 

0 

cos(0 - y) 

0 

-sin(0-y) 

cosH^ 

sin'F w 

0' 

0 

cos <p 

sin 0 

0 

1 

0 

“Siny w 

cos 4^ 

0 

0 

-sin <f> 

cos <p 

sin(0-y) 

0 

cos(0 - y ) 

0 

0 

1 


( 20 ) 


At each time step then the turbulence velocities are 
rotated into body axes with the following equation. 



turbulence 


( 21 ) 


Figure 7 shows the effect of the transformation in the 
case of a complete roll about the longitudinal body axis. 
The roll maneuver begins at the 10 second mark and is 
completed in 24 seconds. The top two traces are the 
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turbulence time histories of the lateral components as 
generated in the turbulence axis system. The bottom 
two traces are the corresponding plots in body axes. 
Notice that the corresponding v and w plots are identical 
for the first and last 10 seconds in the traces. Between 
10 and 34 seconds the roll occurs. Notice at the half 
way point the two v plots and two w plots are equal in 
magnitude but opposite in sign, as they should be for 
180 degrees of roll. 

Turbulence Segments 

The model as described by Tank was developed for 
simulations with the simulation time duration defined a 
priori. The low frequency wave number was then 
computed as the quotient of the time duration and the 
flight velocity. This manner of implementation allows 
one set of Fourier coefficients to generate the entire 
turbulence field. 

When the simulation time duration is not known, the 
low frequency wavenumber must be specified. If the 
duration of the simulation is long enough for the 
vehicle to travel a distance greater then the reciprocal of 
the selected low frequency wavenumber, hereafter called 
the segment distance, the turbulence field would begin 
to repeat. To avoid the repeating pattern, a new set of 
Fourier coefficients must be calculated each time a 
segment ends. 

This changing of Fourier coefficients allows the likely 
possibility that a significant discontinuity could exist 
between segments. To minimize the discontinuity, an 
exponential decay from the last value of the previous 
segment is faired into the following segment according 
to, 

u(t) = u* (t) + ( Au) exp^ At : j (22) 

where u’(t) is the raw turbulence perturbation as 
calculated, Au is the magnitude of the discontinuity 
between the segments, f 2 is the upper frequency limit, 
and At in this instance is the time since the 
discontinuity. The result of equation 23 is that the 
discontinuous jump is faired in over a few seconds in 
most cases. Figure 8 shows a turbulence time history 
sequence in its raw and faired form. In this example the 
minimum frequency of interest was 0.05 Hz, therefore a 
discontinuity occurs every 20 seconds. In Figure 8 they 
are visible in the raw sequence (dotted line) at the 260 
and 280 second points. Notice that the faired sequence 
(solid line) does not follow the discontinuity, but 
quickly merges with the raw sequence shortly thereafter. 


Summary 

In this report Dryden and Kolmogorov based turbulence 
models have been reviewed with attention paid to the 
parameterization of the models. The specific 
Kolmogorov model proposed by Tank has been 
implemented as a series of MATLAB scripts. Features 
that add functionality to the Tank model have been 
developed, implemented and described. A variation of 
the eddy dissipation rate with altitude based on Dryden 
data has been calculated. The eddy dissipation rate 
database has been mapped onto a regular latitude and 
longitude grid and over-land values filled in. Methods 
to account for flight at large vehicle attitude angles have 
been implemented. Finally, a procedure for 
transitioning smoothly across turbulence segments has 
been developed. 
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Table 1. Dryden horizontal and vertical turbulence parameters as functions of altitude and severity. 



Light Turbulence 

Moderate Turb. 

Severe Turbulence 


Altitude 

Horizontal 

Vertical 

Horizontal 

Vertical 

Horizontal 

Vertical 

Horizontal 

Vertical 


<j h 

^v,w 


^v,w 


<L.w 


^V.W 

km 

m/s 

m/s 

m/s 

m/s 

m/s 

m/s 

km 

km 

1 

0.17 

0.14 

1.65 

1.36 

5.70 

4.67 

0.832 

0.624 

2 

0.17 

0.14 

1.65 

1.43 

5.80 

4.75 

0.902 

0.831 

4 

0.20 

0.17 

2.04 

1.68 

6.24 

5.13 

1.04 

0.972 

6 

0.21 

0.17 

2.13 

1.69 

7.16 

5.69 

1.04 

1.01 

8 

0.22 

0.17 

2.15 

1.69 

7.59 

5.98 

1.04 

0.98 

10 

0.22 

0.17 

2.23 

1.73 

7.72 

6.00 

1.23 

1.10 

12 

0.25 

0.18 

2.47 

1.79 

7.89 

5.71 

1.80 

1.54 

14 

0.26 

0.19 

2.62 

1.91 

6.93 

5.05 

2.82 

2.12 

16 

0.24 

0.21 

2.44 

2.10 

5.00 

4.31 

3.40 

2.6 

18 

0.22 

0.21 

2.21 

2.07 

4.07 

3.81 

5.00 

3.34 

20 

0.23 

0.20 

2.26 

1.99 

3.85 

3.38 

8.64 

4.41 

25 

0.27 

0.21 

2.71 

2.09 

4.34 

3.34 

12 

6.56 

30 

0.37 

0.24 

3.73 

2.39 

5.6 

3.59 

28.6 

8.88 


Table 2. Epsilon values derived from Table 1 data. 


Altitude, 

km 

Light Turbulence 
e, m 2 /s 3 

Moderate Turbulence 
e, m 2 /s 3 

Severe Turbulence 
e, m 2 /s 3 

Normalized 
variation, e n 

i 

5.53E-05 

0.0506 

2.09 

2.77 

2 

5.10E-05 

0.0467 

2.03 

2.56 

4 

7.21E-05 

0.0765 

2.19 

3.61 

6 

8.34E-05 

0.0871 

3.31 

4.18 

8 

9.59E-05 

0.0895 

3.94 

4.81 

10 

8.1 IE-05 

0.0845 

3.50 

4.07 

12 

8.13E-05 

0.0784 

2.56 

4.08 

14 

5.84E-05 

0.0598 

1.11 

2.93 

16 

3.81E-05 

0.0400 

0.34 

1.91 

18 


0.0202 

0.13 

1.00 

20 

1.32E-05 

0.0125 

0.06 

0.66 

25 

1.54E-05 

0.0155 

0.06 

0.77 

30 

1.66E-05 

0.0170 

0.06 

0.83 
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Table 3. e at regular latitude and longitude intervals (for proper scaling multiply values by 10' 5 mV 3 ). 
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Figure 1 . Wind and Potential Temperature Spectra as reported by Nastrom 
and Gage (1981). Note: for clarity, the meridional wind and temperature 
spectra have been shifted one and two decades to the right, respectively. 
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Figure 6.- Temperature histories with and without s scaling. 



Figure 1 - Comparison of transverse turbulence in turbulence and body axes 
during 360 degree roll. 
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Figure 8.- Longitudinal turbulence time history showing the effect of ‘fairing’ 
between segments. 
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